Absence of the liquid phase when the attraction is not pairwise additive 
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Abstract 

Recent work on charged colloidal suspensions with very low levels of added salt has suggested that although pairs 
of the colloidal particles repel, clusters of the particles attract. Motivated by this, we study simple model particles 
which have many-body attractions. These attractions are generic many-body attractions and are not calculated for 
any specific colloidal suspension. We find that many-body attractions can stabilise solid phases at low pressures but 
that the liquid phase is either completely absent from the equilibrium phase diagram or present only within a small 
region of parameter space. 



1 Introduction 

Although the behaviour of suspensions of highly charged 
colloidal particles at very low added salt concentrations 
is highly controversial, accurate measurements on isolated 
pairs of the particles shows nothing but a pure repulsion 
yet the particles form crystallites which appear to be 
metastable at close to zero osmotic pressure . This is puz- 
zling, particles interacting via a pairwise additive, purely 
repulsive interaction should only form a solid phase at high 
osmotic pressure, the spheres need to be pushed together 
against their repulsion. The solid phase of hard spheres is 
certainly not metastable at low pressure. Thus the obser- 
vation of the repulsion of isolated pairs, and of the clus- 
tering of larger numbers of spheres is inconsistent with a 
description of the interaction of the colloidal particles via 
a pairwise additive potential of mean force. If we are to 
continue to describe these systems as particles interact- 
ing via some potential of mean force, then we must re- 
lax our restriction to a pairwise additive potential of mean 
force and consider a many body potential We do this 
here. We consider simple, rather generic, many-body at- 
tractions and calculate the phase behaviour they lead too. 
We find that if pairs of particles repel but triplets attract 
that the phase behaviour is qualitatively similar to that 
found for a conventional pairwise additive attraction ex- 
cept that the region within which the liquid phase is found 
is much smaller. If pairs and triplets repel but four particles 
attract then there is no liquid phase at all; the attraction 
merely broadens the coexistence region between the fluid 
and solid phases. The experiments find metastable solid 



phases but no liquid phases, except for one, controversial, 
resuh g. 

Here, we do not address the question of what is the ori- 
gin of a purely repulsive interaction between pairs of parti- 
cles but an attraction between larger numbers of particles. 
We do not attempt to derive the many-body attraction 
from the electrostatic interactions between the charged col- 
loids and the counter and colons. However, there has been 
some work in which the existence, or non-existence, and 
origin of many-body interactions in suspensions of highly 
charged colloids has been considered [|[^. Also, van Roij, 
Dijkstra and Hansen 0| have found both fluid-solid and 
fluid-fluid coexistence within an approximate theory. We 
are simply interested in the consequences, particularly the 
phase behaviour, of many-body attractions. Because we 
have not derived our interaction potentials from the un- 
derlying electrostatic interactions we cannot say anything 
quantitative about the suspensions of highly charged col- 
loidal particles. However, our results show that rather gen- 
erally if pairs repel then a many-body attraction can still 
cause condensed phases to appear but that these condensed 
phases are much more likely to be solid rather than liquid 
phases. This conclusion only relies on the attractions be- 
ing many-body, it does not rely on, for example, whether 
or not the particles are charged. 

In order to keep our model of non-pairwise additive at- 
tractions as simple as possible we simply generalise a simple 
van der Waals-type model. This model is that of particles 
interacting via a hard-sphere repulsion and a long-range 
pair attraction. It has been extensively studied, see Refs. 
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|p|-p^ . We simply generalise the long-range pairwise addi- 
tive attraction to a long-range attraction between n par- 
ticles, i.e., an n-body attraction. This is done in section 
2. In section 3 we apply the standard perturbation the- 
ory to obtain the contribution of the attractions to the free 
energy. Then in section 4 we show and discuss the phase 
behaviour before ending with a conclusion, section 5. 

2 Potential 

Our potential consists of a hard-core interaction which is 
pairwise additive — the hard-sphere potential — plus an 
rt-body attraction. The energy of N particles is 



(1) 



where denotes the N position vectors of the N spherical 
particles, and e„ is a measure of the strength of the n-body 
attraction. It is the sum of two terms. Vhs is the potential 
energy of N hard spheres, given by 

VHs{r'') = l E (2) 

i=l,N j = l.,N 

where 4>hs is the hard-sphere pair potential. 



4>hs{r) 



oo 




r < a 
r > a ' 



(3) 



cr is the hard-sphere diameter. Va is the energy of attrac- 
tion coming from an rt-body attractive potential (f)a, 



K(r^) 



3 E E 0a({|r,-r..|n. (4) 



1=1, N 



The n-fold sum is over the TV!/ (TV — n)l sets of n out of the 
N particles formed by choosing one of the N particles to 
be particle i, then one of the remaining — 1 particles to 
be particle j, ... and then one of the iV — n -I- 1 particles as 
particle /. The dash on the sum over I indicates that the 
sum excludes terms in which I — This overcounts 

by a factor of n! as it counts each interaction nl times, each 
time with the n particles in a different order. For exam- 
ple, for n = 2 it counts the interaction between particles 
numbers 11 and 103, for example, twice, once as i — 11, 
j — 103 and once a.s i — 103, j = 11. For general n the 
number of ways a set of n particles can be assigned to the 
n indices i, j, . . I is n\. The n-body attractive potential 
(/)a is a non-positive function of the set of n(n — l)/2 pair 
separations of the n-particles. This set of scalar separations 



is denoted by {\rj — r^l}". The potential (pa is symmetric 
with respect to exchanging any pair of the n particles. We 
write it as the product of a strength of attraction e„, and a 
function C({|rj — rfe|}") which determines the dependence 
of the energy on the particle coordinates 



MU^j - rfelD = -enC({|rj - rfe|}"), 



(5) 



where C is a non-negative function of the set of pair separa- 
tions of the n particles which we need not specify explicitly 
but which must be long ranged, i.e, it decays to zero over 
some characteristic range much larger than the hard-sphere 
diameter a. For a configuration in which all n(n — 1)/2 sep- 
arations \rj — Tfcl are not much more than the hard-sphere 
diameter, C({ki ~ ^ 1- Then if one or more of the 

separations is increased C tends to zero over a range which 
is much larger than the hard-sphere diameter. The integral 
of C over the positions of all n-particles is proportional to 
the volume V in the thermodynamic limit 



ci{\v,-v,\r)dv- = viy„ 



(6) 



h'n is finite and a constant; it has dimensions of length to 
the power of 3 (n — 1). 

3 Thermodynamic functions 

As the attraction is assumed to be long-ranged, i.e., with 
a range much larger than the hard-sphere diameter cr, we 
may use the van der Waals approximation of approximating 
the free energy by the free energy of a system interacting 
only via the repulsive part of the potential plus the energy 
of interaction due to the attractive part of the potential. 
Here the system interacting only via the repulsions is that 
of hard spheres, whose free energy in the fluid and solid 
phases are accurately given by the approximations of Car- 
nahan and Starhng |13|, and of Hall Q, respectively. So, 



we approximate the Helmholtz free energy by that of hard 
spheres Ahs, plus the energy (not the free energy) evalu- 
ated by assuming that the configurations of the particles 
are unaffected by the attraction, U 

A{N,V,T) = Ahs{N,V,T) + UiN,V). (7) 

Ahs is given by 

A,,s^ -kTlnZhsA^ (8) 

where 

Zhs^^J exp [-msir")] dr^ (9) 
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is the configuration integral for N hard spheres in a volume 
V. A is the thermal volume of a particle. U is given by 

^^_ /K(r^)exp[-/3y,,(r^)] dr^ ^^^^ 



If we substitute our expression for Va, Eq. (Q), into Eq. 
( |lO|) we note that it is the sum of Nl/{N — n)l equivalent 
terms [ITsl. So, 



U 



(N -ny.n\ 



/.^,({|r,--rfe|nexp [-/31^..(r^)] dr^ 

but the n-particle density of hard spheres, Pj^J , and their 
n-particle distribution function, gj^J , are defined as 

\i=l,n I 

N\ /exp [-/?y,,,(r^)] dr^-" 



{N ~ny. 



(12) 



The 1-particle density, p'^j^^ (r) , is not assumed to be uniform 
so that the theory applies to solid as well as fluid phases. 
Using, Eq. ( [l^ ) in Eq. (|ll|), we obtain 



U 



l^ Jc^a ({|r, - r.lD p«(r,) j 5l':Hr")dr". 

(13) 

We now use the long-range of the potential 4>a to simplify 
Eq. (|l3|). As (pa decays to zero only when the separations 
of the particles are much larger than cr, the integral for U 
is dominated by configurations when the spheres are far 
apart, i.e., when all n{n — l)/2 pair separations are much 
larger than a. In a fluid the one particle density is a con- 
stant, p^'^'' = p, and at separations large with respect to a 
the distribution function is close to one, gj^"^ 2± 1. Thus 
in the fluid phase the integrand of Eq. ( [l3| ) is simply ap- 
proximated by p'^4'a- In a solid phase, although there are 
long-range correlations in pj^"'' the one particle density, p')^l , 
averaged over a unit cell is just p. The attractive interac- 
tion between particles has a range much larger than than 
the lattice constant of the lattice and so pa varies little 

across a unit cell and we can regard p\i](jc) as approxi- 
mately constant at its average value, p. Similarly, as we 



change any one of the n position vectors upon which the 
n-body distribution function, p^f^J , depends the density os- 
cillates rapidly over each unit cell but averages to p. So, we 
approximate p^^"'' by p". Then the integrand of Eq. ( |l3| ) 
is the same as in a fluid phase, p^(j)a- So, we approximate 

(") u 
9hs by one, 

solid phases, 



by one, and p^^l] by p in Eq. ( ^ for both fluid and 



(14) 



then using Eqs. d) and (I) we obtain our final expression 
for the energy 



U ^ _!lPl^p^ = -Na,,p:-^kT, (15) 



where pr = po^ is a reduced density and 



an 



(16) 



is inversely proportional to the temperature. 

Using, our expression for the energy, Eq. ( p^ in the van 
der Waals approximation for the Helmholtz free energy, Eq. 
(0) yields 



f3a{p, T) = Patsip) - a„p" 



(17) 



where a = A/N and Uhs = A^s/N. The pressure times a^, 
p, is easily obtained by differentiating the free energy, Eq. 

(Hi 



l3p{p,T) = Pphsip) - {n- l)a„p" 



(18) 



Phs is the pressure of hard spheres times . The chemical 
potential p is then obtained from p = a + p/pr- Knowl- 
edge of the pressure and chemical potential is enough to 
determine the phase diagram. 

For n = 2, Eq. (|l^) reduces to a corrected version of 
the free energy derived by van der Waals 120 years ago 
1^, ^|ll|. By corrected we mean that van der Waals' crude 
approximation for the free energy of particles interacting 
via a strong repulsion has been replaced by the free energy 
of hard spheres. Van Kampen [|l^, and Lebowitz and Pen- 
rose [0 were able to show that if Uhs is exact then Eq. 
( [T^ ) for n = 2 is exact in the limit that the range of the 
attractive potential tends to infinity. For finite range it is 
a reliable and very widely used approximation. 



3 



4 Phase behaviour 

Hard spheres freeze at high densities, i.e., they undergo 
a first-order transition to the sohd phase ||]. The tran- 
sition occurs between a fluid phase at a volume fraction 
77 = {Tr/6)pa^ — 0.49 and a solid phase at the higher volume 
fraction 77 — 0.55 |p^|-p^. This is the only transition in the 
an — limit. The a„ are integrated strengths of attraction 
divided by temperature and so l/a„ is a reduced temper- 
ature. Thus, the a„ = limit is equivalent to the T ^ co 
limit. However, as the temperature decreases a„ increases. 
The first effect this has is to widen the fluid-solid coex- 
istence region due to the fact that the denser solid phase 
has its chemical potential and pressure lowered more than 
the less dense fluid. But at sufficiently low temperature 
the pressure in a single phase can become a non-monotonic 
function of density — the negative term due to attractions 
in the pressure, Eq. ([l8| ) creates a dip in the pressure — 
this is a van der Waals loop and it indicates that there is 
phase separation into dilute and dense phases of the same 
symmetry. The signature of phase coexistence is a van- 
der- Waals loop as our theory is a mean-field theory; it will 
therefore predict incorrect critical exponents. 

For n = 2, we have particles interacting via a hard- 
sphere repulsion and a long-range pairwise additive attrac- 
tion: essentially the model postulated by van der Waals 
in the last century to describe the vapour-liquid transi- 
tion. It, of course, has a vapour-liquid transition, as can 
be seen in Fig. |l| which shows its phase behaviour in the 
density-temperature plane. There is a large temperature 
range over which there are coexisting dilute (vapour) and 
dense (liquid) fluid phases: the temperature at the critical 
point is more than twice its value at the triple point. The 
critical point is the maximum in the vapour-liquid coexis- 
tence curve and occurs at l/a2 — 0.18, and the triple point 
temperature is the temperature at which vapour, liquid and 
solid phases coexist. It marks the lower limit of equilibrium 
vapour-liquid coexistence and occurs at l/a2 — 0.084. 

For n — 3, although there is still equilibrium vapour- 
liquid coexistence, see Fig. ||, the temperature range over 
which it occurs is much reduced. At the critical and triple 
point temperatures, l/aa = 0.20 and 0.17, respectively. 
The ratio between the two temperatures is only 1.2. The 
density at the critical point is also considerably higher than 
for a pairwise additive attraction. 

For n = 4 there is no equilibrium vapour-liquid coexis- 
tence, see Fig. ^ At equilibrium there is only one phase 
transition: the fluid-solid transition. The vapour-liquid 
transition has not disappeared without trace, however. Our 
mean-field pressure in the fluid phase does develop a van 



der Waals loop at low temperature and we can construct 
a vapour-liquid coexistence curve within the fluid-solid co- 
existence region. This is plotted as a dashed curve in Fig. 
||. Note that the critical density is even higher than for 
a 3-body attraction. The disappearance of the equilibrium 
vapour-liquid coexistence is just what is observed when the 
range of a pairwise additive attraction is reduced so that it 
is only about 10% or less of the hard-sphere diameter, see 
Refs. H-H. 

For n larger than 4, the vapour-liquid transition moves 
deeper into the fluid-solid coexistence region, the density at 
the critical point increases with increasing n. If n is made 
extremely large then a van der Waals loop appears in the 
pressure of the solid phase. By n = 15 there is equilibrium 
coexistence between two solid phases of the same symmetry 
but different densities. This is just as has been found with 
a very short-range attraction, a few % of the hard-sphere 
diameter [ p2| . 

The physical interpretation of the phase behaviour is 
straightforward. The energy per unit volume due to an n- 
body attraction is proportional to the number of clusters 
of n particles where all n particles are close enough to each 
other to be within range of the attraction. It is this de- 
pendence which leads directly to the p" dependence of the 
energy per unit volume. So, the larger n is, the smaller 
the number of clusters of n particles at low density and 
the more rapidly the number of clusters grows as density 
increases. When n is small, n — 2, the energy starts to be- 
come significant at a relatively low density and it decreases 
gently, as density squared. This creates a van der Waals 
loop in the pressure at low density. The slow decrease in 
energy as density increases allows the hard-sphere part of 
the free energy to dominate at densities which are not too 
high. This means that the pressure turns up at densities 
below freezing, so allowing a stable liquid to coexist with 
the vapour. For n = 2 the van der Waals loop is below the 
density at which hard-spheres solidify and we see equilib- 
rium fiuid-ffuid coexistence. However, for n > 2, the energy 
only becomes significant at higher densities and it decreases 
rapidly with increasing density. The more rapid decrease 
with density of many-body attractions directly favours the 
solid phase simply because it is more dense. This then 
broadens the fluid-solid coexistence region to such an extent 
that the fluid-fluid transition becomes metastable. Also, as 
the energy only decreases at relatively high density for large 
n, the van der Waals loop is shifted to higher density; ulti- 
mately for n = 6 the critical point is at such high density, 
a volume fraction rj = 0.51, that it is above the density at 
which the particles freeze. It is therefore not observable. 

To summarise: a combination of the increased stabil- 
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ity of the denser solid phase with respect to the less dense 
fluid phase, and the increasing density of critical point drive 
the fluid-fluid transition metastable, then, as we continue 
to increase n, the critical density becomes so large that 
it exceeds the density at which hard spheres crystallise. 
Beyond this point it is not possible to observe fluid-fluid 
coexistence. 



5 Conclusion 

We have shown that many-body attractions have a much 
weaker tendency to induce separation of the fluid phase 
into dilute (vapour) and dense (liquid) phases than do pair- 
wise additive attractions. Indeed when not only pairs but 
triplets of particles interact via a pure repulsion and with 
only four or more particles is there attraction, we found 
no equilibrium coexistence between dilute and dense fluid 
phases. Therefore, we expect that the findings for highly 
charged colloidal particles are general in the following sense: 
if pairs of particles repel then some systems may still have 
an attraction between larger numbers of particles. This 
may cause the fluid-solid transition to broaden greatly and 
produce a solid phase at much lower (osmotic) pressures 
than in a purely repulsive system, however, there will ei- 
ther be only a small range of parameters over which there is 
equilibrium fluid-fluid coexistence or no equilibrium fluid- 
fluid coexistence at all. Of course there are an infinity 
of possible n-body attractions and they will lead to dif- 
ferent behaviour in the same way that different forms of 
pairwise attraction leads to different behaviour. However, 
our results should apply to potentials which are reasonably 
smooth and long-ranged, i.e., with a range of roughly the 
hard-sphere diameter or more. We are also assuming that 
the particles are spherical. 

Finally, there has been a great deal of interest in the 
disappearance of the liquid phase as the potential is var- 
ied; attention has focused on particles in which the range of 
a pairwise additive attraction is very short |l9|-p4[ . But for 
other examples of liquid phases disappearing see Ref . |^ . 
The behaviour we have found as we went from 2 to 3 to 
4-body attractions is qualitatively exactly the same as has 
been found as the range of a pairwise additive attraction is 
made very short. 

It is a pleasure to acknowledge discussions with D. 
Frenkel, B.-Y. Ha and R. van Roij. 
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Figure 1: The phase diagram of hard spheres plus a long-range pairwise additive attraction, n = 2. The a;-axis is a 
reduced density, the volume fraction 77, and the y-axis is a reduced temperature, l/a2- The thick solid curves separate 
the one and two-phase regions. The letters V, L and S denote the regions of the phase space occupied by the vapour, 
liquid and solid phases, respectively. The horizontal thin lines are tie lines connecting coexisting densities. 




Figure 2: The phase diagram of hard spheres plus a long-range three-body attraction, n = 3. See caption of Fig. |^ for 
details. 
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Figure 3: The phase diagram of hard spheres plus a long-range four-body attraction. The letters F and S denote 
the regions of the phase space occupied by the fluid and solid phases. The dashed curve is the coexistence curve for a 
vapour-liquid transition within the fluid-solid coexistence region. See caption of Fig. m for further details. 
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